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HIGHLIGHTS 


•  Novel  approach  to  analytical  modelling  pyrolysis  kinetics  of  oil  shale  is  presented. 

•  The  methodology  is  generic  which  can  model  many  types  of  reaction  kinetics. 

•  The  model  is  efficient,  accurate,  easily  implemented  and  application-independent. 

•  Comparison  with  measurements  and  numerical  simulations  shows  excellent  agreement. 

•  Favorable  generalizations  and  extensions  of  the  modelling  approach  are  presented. 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  30  July  2013 

Received  in  revised  form  2  June  2014 

Accepted  3  June  2014 

Available  online  17  June  2014 


Keywords: 
Analytical  model 
Pyrolysis  kinetics 
Arrhenius  equation 
Oil  shale 


This  paper  presents  a  new  approach  to  analytical  modelling  pyrolysis  kinetics  of  oil  shale  based  on  Arrhe¬ 
nius  equation  which  is  inherently  nonlinear  allowing  no  exact  analytical  solution.  The  proposed  method 
introduces  an  approximation  which  leads  to  a  significant  simplification  of  the  closed-form  solution  and 
the  calculation  and  allows  for  rapid  estimates  of  kinetic  parameters.  The  analytical  model  is  successfully 
validated  by  both  experimental  data  and  the  numerical  model  ( R 2  >  0.99).  In  addition,  the  developed 
methodology  and  model  are  sufficiently  general  which  can  be  applied  to  other  types  of  fuel  resources 
and  chemical  reactions.  To  demonstrate,  some  generalizations  and  extensions  of  representative  kinetic 
models  in  wide-area  applications  are  presented.  The  proposed  modelling  approach  improves  over  the 
state-of-the-art  methods  for  pyrolysis  kinetics  models  in  its  efficiency,  simplicity,  generality  and  applica¬ 
tion-independence.  The  model  is  easy  to  implement  and  computationally  attractive. 

©  2014  Published  by  Elsevier  Ltd. 


1.  Introduction 

As  an  energy  source,  oil  shale  pyrolysis  has  been  produced  over 
hundreds  years  through  heating  the  oil  shale  to  a  temperature  that 
the  kerogen,  the  organic  portion  of  the  oil  shale,  decomposes  into 
gas,  oil  and  coke  [1,2].  Broadly,  two  basic  oil  shale  retorting  pro¬ 
cesses,  in  situ  (underground)  and  ex  situ  (aboveground),  are  used. 
Both  techniques  require  energy  to  heat  the  oil  shale  to  pyrolysis 
temperature.  Modelling  the  pyrolysis  kinetics  of  oil  shale  for  an 
understanding  of  the  decomposition  mechanisms  and  kinetic 
parameters  are  significant  for  designing  and  operating  an  energy 
efficient  pyrolysis  process  because  pyrolysis  is  a  very  complicated 
thermochemical  conversion  process  involving  extremely  complex 
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reactions  and  the  end  products  depend  on  numerous  factors  [3]. 
Hence,  despite  extensive  knowledge  of  the  process,  an  exact  mech¬ 
anism  and  kinetic  modelling  for  oil  shale  pyrolysis  is  difficult  or 
practically  impossible.  Therefore,  many  comprehensive  pyrolysis 
models  are  based  on  Arrhenius-type  kinetic  equation  for  interpret¬ 
ing  fundamental  mechanism  of  the  entire  pathway  or  specific 
mechanisms  within  the  pathway.  Arrhenius  model  has  a  long  his¬ 
tory  of  use  as  pyrolysis  kinetics  and  revisiting  the  mechanisms  for 
pyrolysis  and  the  relevant  models  so  far  developed  demonstrate  its 
generality  which  has  proven  to  be  a  valuable  analysis  tool  with 
wide  applications  [3-8  . 

Arrhenius  kinetics  for  decomposition  is  modelled  based  on  the 
two  most  important  and  frequently  used  parameters:  activation 
energy  and  pre-exponential  factor  (e.g.  [4  ).  Since  an  accurate  Arrhe¬ 
nius  model  equation  is  inherently  nonlinear  in  nature  allowing  no 
exact  analytical  solution,  numerical  techniques  become  necessary 
to  find  the  solutions  and  the  parameters  (e.g.  [9  ).  The  underlying 
mechanisms  remain  elusive  with  numerical  approaches. 
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The  analytical  model,  on  the  other  hand,  has  an  advantage  that 
it  can  provide  a  clear  description  of  the  primary  kinetics  and 
parameters  for  the  underlying  mechanisms.  It  is  suited  to  under¬ 
stand  the  outcomes  and  to  interpret  the  experiments.  Extensive 
analyses,  therefore,  have  been  performed  to  explore  thermal 
decomposition  kinetics  and  to  determine  kinetic  parameters  for 
pyrolysis  kinetics  of  oil  shale  as  well  as  for  biomass  using  analytical 
models  (see,  for  example,  [10,11],  for  general  reviews).  Two  basic 
principles  regarding  the  methodology  can  be  identified  among 
the  analytical  models:  integral  and  differential  methods  based  on 
Coats-Redfern  12]  and  Friedman  13]  approaches  and  their  mod¬ 
ifications  [6,14-16].  Typically,  Thermogravimetry  (TG)  and  Differ¬ 
ential  Thermal  Analysis  (DTA)  experiments  are  conducted  under 
either  isothermal  or  non-isothermal  condition.  The  TG/DTG  data 
are  analyzed  and  the  kinetic  parameters  are  adjusted  based  on 
the  kinetic  model  equation  using  goodness-of-fit  assessment 
[17].  Graphical  methods  are  often  used  to  evaluate  goodness-of-fit. 

Although  simple,  an  obvious  disadvantage  with  these  methods 
is  that  it  describes  observations  of  the  data  only  without  showing 
the  underlying  mechanism  by  which  the  data  are  produced. 
Another  relevant  disadvantage  which  is  directly  tied  to  this  is  that 
graphical  methods  generally  suffer  from  poor  statistical  properties 
(see,  for  example,  [10  ).  Attempts  have  been  made  to  improve 
these  methods  and  majority  of  the  improvements  are  made  to  bet¬ 
ter  approximate  the  Arrhenius  temperature  integral  (see 
[4,6,18,19]  for  a  general  review).  The  temperature  integral  is 
always  encountered  which  permits  no  analytical  solution.  The 
approximation  of  the  temperature  integral  is  usually  done  through 
predefined  formula,  which  is  undoubtedly  related  to  the  unknown 
kinetic  mechanisms.  Accuracy  is  often  sensitive  to  the  form  of  the 
assumed  kinetic  mechanisms  [6,20,21  ,  which  presents  a  major 
challenge  in  kinetics  analysis.  This  requires  prior  knowledge  of 
the  kinetic  mechanisms  which  is  often  inaccessible  in  practice. 

This  paper  presents  a  new  analytical  methodology  to  address 
such  an  issue  to  improve  upon  past  analytical  models.  We 
approach  the  temperature  integral  problem  differently  than  exist¬ 
ing  models  in  that  the  proposed  model  does  not  make  any  assump¬ 
tions  on  the  temperature  integrals  or  kinetic  mechanisms.  The 
developed  model  provides  a  closed-form  expression  of  the  decom¬ 
position  kinetics  based  on  kinetic  parameters.  Experimental  and 
numerical  validations  show  an  exceedingly  high  accuracy  of  the 
model  (R2  >  0.99). 

One  significant  aspect  of  the  developed  modelling  method  is  the 
straightforward  formalism  for  the  final  model  solution  which  rep¬ 
resents  a  significant  simplification  of  the  closed-form  solution  and 
calculation  and  allows  for  rapid  estimates  of  kinetic  parameters. 
The  proposed  approach  is  highly  efficient,  simple,  accurate,  yet 
universal  which  can  be  easily  modified  and  extended  to  accommo¬ 
date  other  reaction  kinetic  features.  To  demonstrate,  we  present  its 
generality  and  practicality  through  representative  wide-area  appli¬ 
cations  and  show  that  many  popular  existing  models  can  be 
expressed  as  instances  of  our  model.  Therefore,  although  the  paper 
concentrates  on  oil  shale  pyrolysis,  the  proposed  methodology  will 
be  of  great  interest  to  researchers  and  engineers  for  modelling 
other  types  of  chemical  reactions  for  fuels  [22  .  Finally,  the  imple¬ 
mentation  is  extremely  simple  which  can  be  achieved  with  a  few 
lines  of  code  or  with  spreadsheets  or  even  with  a  pencil  and  paper. 

2.  Methods 

2.2.  Kinetic  model 

The  kinetics  of  kerogen  decomposition  to  shale  oil  is  described 
in  two  steps  [23,24]: 

kerogen  ->  bitumen  +  gas  ->  oil  +  gas  +  residue 


and  modelled  as  an  n-order  reaction  equation  [25]  for  retorting 
process 

dX 

Tt=k(\-Xf  (1) 

where  X  is  the  weight  loss  or  mass  fraction  or  conversion  (%) 
defined  as 


w0  -  w t 

W0  -  Wf 

w0:  initial  weight; 
wt:  weight  at  time  t; 
wf.  final  weight; 

k:  rate  coefficient  given  by  Arrhenius  equation: 


k  =  A  exp  (T  T j  (3) 

A:  pre-exponential  factor; 

E :  activation  energy; 

R:  gas  constant; 

T:  temperature. 


Assuming  a  constant  heating  rate  p  =  an  integration  of  Eq.  (3) 

results  in  a  plot  of  In  (-  ln(|2~X))  or  (in  ")  if  n  >  l)  against 

1  IT  being  a  straight  line.  The  kinetic  parameters  are  determined 
using  its  slope  and  intercept.  This  graphical  method  is  based  on 
Coats-Redfern  method  which  is  the  most  commonly  used  param¬ 
eter  estimation  approach.  Despite  of  ease  of  use,  graphical  method 
has  some  limitations  as  described  in  Section  1.  In  the  following  we 
propose  a  new  model  to  overcome  these  limitations. 

2.2.  The  proposed  approximate  analytical  model 


Assume  n  =  1.  As  we’ll  see  later  that  this  assumption  is  for  a 
simple  illustration  only  and  the  reaction  order  is  not  limited  to  lin¬ 
ear  one.  Since  a  constant  heating  rate  p  =  g  is  the  most  common 
case  in  applications,  Eq.  (1)  leads  to 


dX 

1  -X 


Integration  of  Eq.  (4)  gives 


dX 
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The  right-hand  integral  cannot  be  calculated  analytically.  We 
propose  its  approximation  as  follows: 
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RT 

T  2 RT 


dT  = 


RT 


exp  - 


RT 


exp 


dT 


RT 


(6) 


Rearrange  of  the  terms  containing  integrals  together  leads  to 
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*T  2 RT 

1  +  -^ 


exp  - 


RT 


dT 


RT 


exp  - 


RT 


(7) 


For  a  small  temperature  change  AT,  1  +  is  approximately 
constant.  Therefore,  Eq.  (7)  can  be  approximated  as 


RT( 
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RTr 


1  + 


2RTn 
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and  Eq.  (5)  becomes 
-  ln(l  —  X)  =  —  ln(l  -  X0)  +  yT  —  yTo 

..  ^^exp(-i)  (9) 

Ir-p  1+2|i 

In  summary,  the  weight  loss  X  can  be  computed  based  on  the 
following  simple  closed-form  formula: 

ln(l  -Xi)  =  -yh  +yTo 

ln(l-X,)  =  ln(l-X,_1)-7T|+7TM,  i  =  1,2,...  (10) 

It  can  be  seen  that  the  proposed  model  is  a  systematic  approach 
which  is  simple,  efficient  and  fully  automatic.  Extensions  of  it  will 
be  presented  later. 

3.  Results  and  discussion 


and  showed  good  reproducibility.  For  the  purpose  of  validation  and 
verification  of  the  developed  analytical  model,  we  will  focus  on  the 
TG/DTA  measures  under  an  air  atmosphere.  The  detailed  descrip¬ 
tion  of  the  experiment  was  provided  in  [26  .  Fig.  1  shows  the 
measurements. 

Fig.  2  shows  the  comparison  of  the  model  predictions  with  the 
measurements.  The  difference  in  temperature  between  two  time 
periods  was  taken  as  AT  =  25  °C.  The  model  accuracy  assessment 
is  presented  in  "able  1,  where  model  comparison  with  measure¬ 
ment  and  numerical  simulation  data  on  the  basis  of  MAE,  RMSE, 
MSE  and  R 2  are  reported.  These  error  metrics  test  the  differences 
between  the  proposed  model  and  the  numerical  model  as  well  as 
how  well  the  developed  model  fits  the  measurement.  Zero  value 
of  MAE,  RMSE,  MSE  and  R 2  value  of  +1  indicate  a  perfect  fit.  The 
predicted  weight  losses  were  in  excellent  agreement  with  mea¬ 
surements  across  all  these  comparison  metrics  (R2  >  0.99). 


This  section  presents  the  validation  of  the  developed  analytical 
model  against  both  the  measurement  and  the  numerical  simula¬ 
tion.  The  following  frequently  used  performance  measures  are 
employed  to  assess  the  model  accuracy:  mean  absolute  error 
(MAE),  root  mean  squared  error  (RMSE),  mean  squared  error 
(MSE)  and  coefficient  of  determination  (R2)  defined  as 


MAE  —  ^  ^  ^j—\  ^  measured  /simulated  Y predicted 


RMSE  = 


1 


N^i=  1  (y measured /simulated  Y; predicted ) 


1  1/2 


MSE  —  ^  ^._1  (V measured / simulated  Y predicted) 


(11) 


3.2.  Model  validation  against  numerical  simulation 

To  further  test  the  accuracy  of  the  proposed  model,  a  MATLAB 
program  was  developed  to  simulate  the  exact  model  Eq.  (1)  with 
£  =  60kjmor1  and  A  =  34s-1.  Different  temperature  intervals 
(AT  =  2  °C,  30  °C,  50  °C)  and  0-800  °C  temperature  range  were  used 
for  generating  the  temperature  versus  weight  loss  data.  Fig.  3  dis¬ 
plays  the  comparison  results. 

Excellent  agreement  was  obtained  between  the  model  and  the 
simulation.  For  temperature  interval  AT  =2  °C,  two  curves  are 
nearly  identical.  Very  good  accuracy  was  obtained  even  when  the 
temperature  interval  AT  was  chosen  as  50  °C.  These  results  vali¬ 
date  the  applicability  of  the  proposed  model.  Additionally,  esti¬ 
mated  activation  energy  E  =  60.12  kj  mol-1  and  pre-exponential 
coefficient  A  =  34.03  s-1  were  obtained  which  are  very  close  to 
their  real  values.  In  the  next  subsection,  we’ll  present  how  the 


N7L i  Y measured /simulated^ predicted  X)/=l  ^ measured /simulated T^i— i  ^ predicted 

2 

2:3 

w 

1 _ 1 

2 

measured  /  simulated 

y)/— i  y measured  /simulated) 

NJ2i= 1  y predicted 

predicted) 

where  the  subscripts  ‘measured/simulated’  and  ‘predicted’  present 
the  measured  or  the  simulated  and  the  model  predicted  quantities 
and  N  the  total  number  of  the  points.  After  the  validation  procedure, 
we’ll  demonstrate  some  extensions  of  the  developed  model  which 
are  considered  as  most  commonly  used  models  in  reaction  kinetics 
analysis. 

3.1.  Model  validation  against  experiment 

Firstly,  experimental  data  were  collected  from  the  measure¬ 
ment  of  oil  shale  samples  from  Huadian,  China.  The  initial  weight 
of  the  samples  was  about  4  mg  (density  1570  kg  m-3).  TG/DTA 
was  conducted  for  the  samples  using  a  Netzsch  STA  449C  thermal 
analyzer  system  (Germany)  under  non-isothermal  conditions  of 
both  ambient  atmosphere  and  nitrogen  atmosphere  (fluxed  by 
nitrogen  with  a  flow  rate  of  50  cm3  min-1).  TG/DTA  provides  the 
precise  measure  of  the  mass  change  in  regards  to  temperature 
change,  which  is  used  in  analyzing  mass  change,  thermal  resisting 
properties  and  reaction  velocity  in  samples.  Before  the  experiment, 
the  TG/DTA  instrument  was  calibrated  and  the  samples  were  finely 
ground  through  a  63  pm  sieve  to  eliminate  heat  and  mass  transfer 
effects  during  pyrolysis.  During  the  experiment,  the  heating  rate 
was  kept  as  10°Cmin-1  starting  from  ambient  temperature  and 
rising  to  800  °C.  All  the  measurements  were  repeated  several  times 


kinetic  parameters  are  determined  using  the  proposed  analytical 
model. 

3.3.  Model  extension  and  discussion 

The  model  can  be  extended  to  accommodate  more  complex  fea¬ 
tures.  We  will  present  several  of  them  that  are  considered  as  espe¬ 
cially  important  and  popular  in  pyrolysis  kinetics  studies  and  show 
that  these  models  can  be  expressed  as  instances  of  our  model. 

3.3.1.  Higher-order  reaction 

For  any  n-order,  n  >  1,  reaction,  the  approximate  analytical  solu¬ 
tion  Eq.  (10)  can  be  modified  (details  omit)  and  weight  loss  X  is 
extended  as  the  following  general  and  simple  formula: 


(1-X!)1-'7 

11  -  1 


-yTl  +  yTo 


(l-X,-)1-"  (1 


11  -  1 


n  - 1 


yTi  +  yTl_ ,,  >  =  i,2, 


with  yT  = 


A  exp(-j^) 

P  l+¥ 


n  >  1 


(12) 
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Temperature  (°C)  Temperature  (°C) 

Fig.  1.  TG/DTG  measures  of  the  oil  shale  sample  (heating  rate  10  °C  min-1). 


Fig.  2.  Comparison  of  model  predictions  and  measurements. 


Table  1 

Model  accuracy  metrics  for  comparison  of  the  proposed  model  with  measurement 
data  and  numerical  simulation  data. 


Comparison  with 

MAE 

RMSE 

MSE 

R2 

Measurement  data 

0.243 

0.280 

0.081 

0.994 

Numerical  simulation  data 

0.001 

0.003 

0.000 

0.998 

Temperature  (°C) 

Fig.  3.  Comparison  of  model  predictions  and  numerical  simulations. 


employed  to  obtain  the  simulated  data  {Tu  Xu}  i  =  0,. .  .,N.  The  error 
function  is 

X(A,E)=  ^(X,0-Xu)2 

1=1. ..N 

The  parameter  estimates  of  A  and  E  can  be  computed  or 
adjusted  iteratively  by  minimizing  the  error  function  through  gen¬ 
eralized  reduced  gradient  method.  We  found  that  the  iteration 
converged  rapidly.  In  most  calculation  cases,  only  a  few  iterative 
steps  were  needed.  Besides,  many  software  programs  are  available 
for  solving  such  problems,  such  as  MATLAB’s  lsqnonlin  function  in 
Optimization  Toolbox. 

To  assess  the  accuracy  of  the  proposed  approach  we  predefined 
kinetic  parameters  and  applied  numerical  model  to  generate  the 
temperature  versus  weight  loss  data.  Kinetic  parameters  were  esti¬ 
mated  with  the  proposed  method  and  the  comparison  results  are 
displayed  in  Table  2.  The  predefined  parameters  are  listed  as  ‘true’ 
values  of  the  parameters  in  Table  2.  High  accuracy  was  obtained. 
The  estimates  are  nearly  the  same  as  the  ‘true’  values. 

3.3.3.  Temperature  dependent  kinetic  parameters 

In  majority  of  studies  the  kinetic  parameters  of  activation 
energy  and  pre-exponential  coefficient  are  considered  to  be  tem¬ 
perature  independent  (see  all  cited  papers).  However,  experiments 
have  suggested  that  temperature  is  an  important  determinant 
parameters  for  pyrolysis  kinetics  and  accurate  estimates  of  their 
temperature  dependence  are  critical  to  understanding  thermal 
activation  processes  [27  .  However,  classic  graphical  methods  can¬ 
not  unambiguously  measure  and  estimate  their  temperature 
dependence  due  to  the  mathematical  impossibility  [27]. 

To  accommodate  temperature  dependence  in  kinetic  parame¬ 
ters,  the  proposed  analytical  model  is  extended.  We  divided  the 
0-800  °C  temperature  range  of  the  experimental  TG/DTG  data  into 
equal  intervals,  each  of  which  contains  20  °C  variation.  Using  Eq. 
(10),  E  and  T  can  be  calculated  for  each  temperature  interval.  The 
parameters  set  E(T)  and  A(T)  was  then  obtained  across  all  intervals 
and  shown  in  Fig.  4.  As  the  temperatures  were  divided  into  suffi¬ 
ciently  small  intervals,  the  discrete  values  of  the  kinetic  parame¬ 
ters  should  be  expressible  as  continuous  functions  of 
temperatures.  Fig.  4  indicates  that  the  variation  of  activation 
energy  with  temperature  for  oil  shale  at  constant  heating  rate  is 
so  small  that  it  can  be  neglected.  We  may  conclude  that  for  con¬ 
stant  heating  rate,  it  is  reasonable  to  approximate  the  activation 


3.3.2.  Kinetic  parameters 

Using  the  analytical  model  Eq.  (10)  or  Eq.  (12)  the  kinetic 
parameters  of  activation  energy  and  pre-exponential  coefficient 
can  be  easily  estimated  analytically  based  on  the  TG/DTG  measure¬ 
ment.  Assume  n  =  1  for  illustration  and  consider  temperature  ver¬ 
sus  weight  loss  data  {Ti,Xi>0}i=o,.  ..>n  from  TG/DTG  measurement. 
Starting  from  initial  values  of  parameters  A\  and  £i,  Eq.  (10)  is 


Table  2 

Comparison  of  the  known  true  activation  energies  and  estimated  activation  energies. 


True  E  (kj  mol  ^ 

Estimated  E  (kj  mol 

60 

60.12 

150 

150.10 

300 

300.22 

186 


X.  Lii  et  al./Fuel  135  (2014)  182-187 


100  r 

90  ■ 

80  - 

o 

S 

70  - 

>—9 

60  - 

OX 

u 

<u 

50  - 

s 

c 

40  - 

o 

30  - 

•  pm 

o 

20  - 

< 

10  - 

0 


0  100  200  300  400  500  600  700  800  900 

Temperature  (°C) 


Fig.  4.  Effect  of  temperature  variation  on  activation  energy  at  constant  heating  rate. 


Fig.  5.  Effect  of  heating  rate  on  weight  loss. 

energy  as  being  independent  of  the  temperature,  as  most  reports 
assume. 

Furthermore,  note  that  this  proposed  parameter  estimate 
method  can  be  also  used  to  estimate  both  the  reaction  order  and 
constant  approximation  of  the  kinetic  parameter  through  calculat¬ 
ing  the  mean  of  the  set  of  estimated  parameter  values. 

3.3.4.  Investigation  on  kinetic  mechanisms 

After  validation,  the  proposed  model  can  be  employed  to  inves¬ 
tigate  various  kinetic  mechanisms  of  pyrolysis,  for  example,  the 
effects  of  kinetic  parameters,  such  as  the  initial  weight,  heating 
rate  and  the  final  heating  temperature.  Fig.  5  shows  such  a  simula¬ 
tion  example.  The  results  demonstrate  that  the  peak  rates  of 
weight  losses  shift  towards  higher  temperatures  as  the  heating 
rates  increase.  However,  heating  rates  have  negligible  effects  on 
the  simulated  final  weight  losses. 


the  number  of  temperature  intervals.  Similar  deduction  to  an  n- 
order  reaction  can  be  obtained  also. 

4.  Discussion 

Since  an  accurate  Arrhenius  type  equation  for  modelling  ther¬ 
mal  decomposition  kinetics  of  pyrolysis  is  inherently  nonlinear 
allowing  no  exact  analytical  solution,  the  proposed  method  intro¬ 
duced  an  approximation  to  the  solution.  The  approximation  is  min¬ 
imal,  accurate,  efficient  and  simple  which  precisely  captures  the 
thermochemical  mechanisms  of  oil  shale  pyrolysis.  This  paper 
makes  three  main  improvements  over  the  state-of-art  pyrolysis 
models:  (1)  the  innovative  methodology  for  modelling  pyrolysis 
kinetics  based  on  analytical  approximation  which  allows  various 
extensions  to  accommodate  more  complex  features  in  reaction 
models;  (2)  the  analytical  model  which  represents  a  significant 
simplification  and  efficiency  of  the  kinetic  calculation  and  analysis; 
(3)  the  robustness  and  high  accuracy  of  the  approach  and  model 
which  capture  the  thermal  decomposition  kinetics  of  pyrolysis 
and  allow  for  analytical  estimates  of  kinetic  parameters.  Experi¬ 
mental  and  numerical  validations  show  high  accuracy  of  the  model 
(. R 2  >  0.99).  Implementation  is  extremely  simple.There  is  one  ques¬ 
tion  regarding  to  the  proposed  model.  The  model  accuracy  depends 
on  the  temperature  interval  AT.  However,  a  relatively  large  tem¬ 
perature  increment,  for  example,  at  least  25  °C  (or  even  50  °C, 
see  Fig.  3),  can  be  chosen  for  accurate  results.  The  maximum  allow¬ 
able  temperature  interval  for  a  desired  accuracy  needs  further 
investigation. 

5.  Conclusions 

This  paper  reports  a  novel  method  for  modelling  pyrolysis 
kinetics  of  oil  shale  which  presents  several  methodological  contri¬ 
butions  to  the  current  literature.  Results  approve  that  the  proposed 
efficient  model  is  suitable  for  simulation  of  pyrolysis  kinetics  of  oil 
shale.  Importantly,  the  new  approach  is  based  on  general  ideas  and 
thus  easily  applicable  to  other  more  complex  kinetic  models.  The 
developed  methodology  may  be  of  assistance  in  modelling  any 
type  of  fuel  resources  and  reaction  kinetics  that  follow  Arrhenius 
equation.  The  developed  model  is  simple,  accurate,  efficient,  appli¬ 
cation-independent  and  computationally  attractive.  There  is  a 
need  for  such  models  in  fuel  and  energy  technology. 
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